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APPLICATION OF A MODIFIED FAST 
FOURIER TRANSFORM TO CALCULATE 
HUMAN OPERATOR DESCRIBING FUNCTIONS 

By Richard S. Shirley 
Electronics Research Center 


SUMMARY 

A version of the fast Fourier transform (FFT) is used in a 
hybrid computer program to permit processing of tracking data to 
yield the human operator's describing function almost immediately 
after the period of data-taking. The use of the FFT allows the 
final calculation time required to process 216 seconds of track- 
ing data to be reduced to 3 seconds from the 10 minutes previous- 
ly required on the same computer. The algorithm used permits the 
bulk of the analysis of the data to be performed while the data 
are being taken, and does not require all the data to be present 
in core before processing begins. 


a 


xk 


TABLE OF SYMBOLS 


the index of summation for the additive portion of 
the FFT 

a constant which weights the sinusoids composing the 
system input, see Table I 

the real part of the truncated Fourier transform of 
x(t) at the frequency given by 



x(t) cos (w^t) dt 


B 


xk 


the imaginary part of the truncated Fourier transform 
of x(t) at the frequency , given by 



sin (w^t) dt 


c a subscript referring to the output of the human operator 

at the control stick (see Figures 1 and 4) 


c(nAt) the data samples taken at the human operator's output 



TABLE OF SYMBOLS (cont'd) 


an integer devisable by 4 used to determine the input 
frequencies 

e a subscript referring to the input to the human operator 

at the oscilloscope (see Figures 1 and 4) 

e(nAt) the data samples taken at the human operator input 

FFT Fast Fourier Transform 

F x (w k ) the truncated Fourier transform of x(t) at the frequency 
to^, given by 



x(t)e i dt 


h a subscript used to denote frequencies between the input 

frequencies, equal to 1, 2, 3,... 

i a subscript referring to the system input (see Figures 1 

and 4) 

i(nAt) the system input equals 
14 

X) A, sin (u.nAt) 
k=l K K 


j the square root of -1. 

k a subscript used to denote the input frequencies, equal 

to 1, 2, 3,... ,14 

m(nAt) the data samples taken at the system output 

n the index of summation for the multiplicative portion of 

the FFT 


N 

T 

V“) 



(to) 


the number of data samples taken, 10,800 

the period of data-taking, equal to 216 seconds 

the dynamics of the controlled element (see Figures 1 
and 4 ) 

the linear portion of the quasi-linear describing function 
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(N/D k ) -1 


6 k 

Y k (D k /4) -1 

At the time increment between interrupts, and hence the 

time between data samples, equals .02 sec 

u) k the frequencies of the sinusoids comprising the system 

input, see Table I 

frequencies between the ' s 

$i c (to) the cross power spectral density between the human 
operator's output and the system input 

$i e (w) the cross power spectral density between the human 
"" operator's input and the system input 

$ nn ( w) the continuous power spectral density of the human 
operator's remnant 

INTRODUCTION 

Only recently have dynamic models of the human operator 
been used effectively in the design of man-vehicle systems. This 
is due partially to a lack of understanding of the human operator 
and also to the difficulty and expense of experimentally deter- 
mining values for the various parameters of existing models. 
Improvements in computers and computational techniques are over- 
coming these difficulties, and already it is possible to bring 
about significant .improvements in a man-vehicle system through 
the use of pilot models in preliminary design (refs. 1 and 2). 
This paper describes a computational technique which reduces 
greatly the cost of obtaining values permitting the use of a 
current pilot model, i.e., the quasi-linear describing function. 

One way to characterize the behavior of a human operator 
in a continuous tracking task is by a quasi-linear describing 
function, which consists of a linear describing function and a 
remnant. The linear describing function is the average frequency 
response of the human operator, i.e., his amplitude ratio and 
phase as a function of frequency. The remnant, characterized by 
a continuous power spectral density, is that portion of the human 
operator's output which is not linearly correlated with his in- 
put. The total output of the human operator is the sum of the 
remnant and the output of the linear describing function* (see 
Figure 1) . 

*Examples of human operator describing functions are shown in 
Figures 2 and 3 . 
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HUMAN OPERATOR 
DESCRIBING FUNCTION 










A direct way to measure describing functions in the labora- 
tory involves the use of a hybrid computer and the method of 
Fourier coefficients. The method of Fourier coefficients has been 
extensively investigated and is described in detail (ref. 3) . It 
will be briefly outlined here for completeness . The human opera- 
tor is placed in a control loop, possibly as shown in Figures 1 
and 4. The system input, a sum of sinusoids of known amplitude, 
phase, and frequency is updated every At seconds; simultaneously, 
data are taken at the human operator's input and output. At the 
end of T seconds, the sampled values of the human operator's in- 
put and output are processed as follows: 


A ck = 

A . N 

c(nAt) 

n=l 

cos 

( w^nAt) 

(1) 

B ck = 

At N 

^£c(nAt) 

n=l 

sin 

(ai^nAt) 

(2) 

A ek = 

M- N 

-tjt 12 e (nAt) 
n=l 

cos 

(to k nAt) 

(3) 

B ek " 

e (nAt) 

1 n=l 

sin 

(oyiAt) 

(4) 

W 

= A ck " jB 

ck 

i 

/ 

f 

(5) 

W 

= A ek " ^ B 

ek 

/ 

/ 

/ 

(6) 


w 





e vu V 


(7) 


/Y (a) ) = -tan" 1 + tan" 1 ^ 

/ P k A_,_ A 


ck 


ek 


( 8 ) 


/ 
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Figure 4.- Flow diagram of the human operator in a compensatory 
tracking system. 
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P k 
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2 2 

A ek + B ek 


1/2 
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$ _(w, ) 

cc h 


2ttT l F c ^b) I 


( 10 ) 


4 > (<0 = 

nn n cc h 


1 + Y Y ( a), ) 
p c h 


( 11 ) 


where the w 's are the input frequencies, and the to, 's lie between 
the w k 's. 

This paper describes how a version of the fast Fourier 
transform (FFT) is used to compute human operator describing 
functions, or more specifically, how a version of the FFT is used 
to solve Eqs . (1) through (4) , while the data samples, c(nAt) and 

e(nAt), are being taken. The FFT is an algorithm which greatly 
reduces the time required to calculate the truncated Fourier 
transform, or periodogram, of a sampled time signal. The savings 
are obtained by replacing calculations which involve trigonometric 
functions or multiplications with simple additions. The replace- 
ment is accomplished by taking advantage of the symmetries of the 
sine and cosine functions, and by further taking advantage of 
relationships between the frequencies at which the Fourier analy- 
sis is performed. 


THE VERSION OF THE FFT USED 

The version of the FFT used takes advantage only of the 
symmetries of the sine and cosine functions. It does not take 
advantage of the relationships among the frequencies at which 
the Fourier analysis is performed. By not using the complete 
version of the FFT, it becomes possible to perform the bulk of 
the data-processing during the At seconds between interrupts 
while the experiment is still in process. The requirement that 
the data be in core before processing, or even that the data fit 
in core, is avoided. The following derivation of the algorithm 
used will make this point clearer. It should be noted that before 
the FFT was used it was not possible to perform the calculations 
between interrupts because the computation time required was over 
two and a half times greater than that which was available. 
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It is desired to evaluate Eqs. (1) through (4) using a 
digital computer. In order to permit the use of the FFT , the 
input frequencies, will be restricted to 


- 2tt 
W k ~ D k At 


where the Dj, are chosen from 4, 8, 12, 16, etc. The method of 
Fourier coefficients further requires that the ratio N/d^ be an 
integer (where N is the number of data samples taken at intervals 
At) . The derivation for A e ^ and is identical to the deriva- 
tion which follows for and B ck . 


Using the identities 



sin 

<N 

+ 

CD 

= sin 0 , and 


cos 

(0 + 2 7T ) 

= cos 0 

or 

sin 

(m k At) = 

sin ^(aD k + 1) w k AtJ 


cos 

(w k At) = 

cos [( aD k + 1) Wj^AtJ 

permits 

Eqs. (1) 

and (2) 

to be rewritten as 


a — 0,1, 2, 3,... 




A , °k 

X 
c a 

i _ 

A ck 

= 

At y 
T 

cos (w k nAt) X) c 


n=l 

a=0 



..A 

e k 

B ck 

= 

At 

T 

sin (u> k nAt) ^2 c 


n=l 

a=0 

(N/D 

k } 

-1. The identities 

sin 

(0 

-it) = 

-sin 0, and 

cos 

(0 

- tt) = 

cos 0 


+ aD. 


+ aD, 


>] 




( 12 ) 


(13) 
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or 


sin 


( M 

a). At 
k 

( n ' -)\ 


p 

/ D, \” 

cos 

o) k At 

> n ""r) 


_ 

\ /j 


= -sin (w^nAt) 


= cos (w^nAt) 


permit Eqs. (12) and (13) to be written as 


ck 


D k /2 

Atr> 

T 2-f 


- c 


cos (oj^nAt) 52 / c^At (n + aD k )J 
a=0 ( 


D k / 2 r 


B 


ck 




n=l 




sin (ci) k nAt) 52 <c[At(n + aD k )] 


- c 


_ At ( n + -T + aD k) 


Finally, the identities sin (-6) = -sin 6, cos (-0) = cos 0, 
sin (tt/2) = cos (tt) = 1 , and sin (tt) = cos (tt/2) = 0 permit 
Eqs. (14) and (15) to be written as 



(14) 


(15) 


(16) 
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where y^ = (D^/4) -1. Equations (16) and (17) represent the 
algorithm used in the hybrid program. The summation over "a" is 
performed between interrupts during the experiment and is called 
the "additive portion" of the FFT. At the end of the data-taking 
period, the summation over n (called the "multiplicative portion" 
of the FFT) and the calculation of the human operator's describing 
function [using Eqs. (8) and (9)], can be performed in less than 
three seconds. 

The hybrid computer program is written in a Fortran IV lan- 
guage which includes hybrid commands. The program is listed in 
Appendix A. Table I lists the values of the experimental param- 
eters, including those which characterize the system input. Fig- 
ure 5 is a flow diagram of the additive portion of the FFT. 

Figure 6 shows a flow diagram of the hybrid program, and lists the 
time taken by each part of the program, both for the FFT version 
and for the version written the old way [directly computing Eqs. 

(1) through (4)]. As shown in Figure 6, the FFT permits a saving 
of nearly ten minutes per run, effectively reducing the run time 
to the time required to take the data and print the results. 


RESULTS 

An initial check of the hybrid program was made by taking 
measurements across known filters. The results shown in Figures 
7 and 8 are quite accurate, and are repeatable. 

Measurements were then taken of the author's tracking per- 
formance in a control loop, as shown in Figures 1 and 4. Ten 
runs were made with each of the controlled elements, 1/s and 
1/s , The describing functions shown in Figures 2 and 3 are com- 
parable with established results (ref. 3). 
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TABLE I 


PARAMETER VALUES USED FOR THE HYBRID PROGRAM 


k 

A k 

volts 

“k 

rad/sec 

n 

A k 

volts 

w k 

rad/sec 

1 

.2 

26.18 

8 

-1. 

1.309 

2 

- .2 

15.71 

9 

1 . 

. 8727 

3 

.2 

8.727 

10 

- 1 . 

.5818 

4 

- .2 

6.545 

11 

1 . 

. 4363 

5 

.2 

4.363 

12 

- 1 . 

.2909 

6 

- 1 . 

2.618 

13 

1 . 

.1745 

7 

1 . 

1.745 

14 

- 1 . 

.1164 


At = time between interrupts = .02 sec 


T^ = warm-up time before data-taking = 24 sec 
T = period of data-taking = 216 sec 

14 

i (nAt) = system input = V A, sin (co.nAt) 

k= 1 k K 


No comparison is made between results for the programs with 
and without the FFT (on Figures 2, 3, 7, 8) because the results 
are identical, as is shown analytically in the derivation of 
Eqs . (16) and (17). The comparison between the computation times 

for the two programs (Figure 6) however, indicates the substan- 
tial savings obtained by using the FFT. The only penalty paid 
for the reduced computational time is an increase in the com- 
plexity of the written Fortran program, as shown in Appendix A. 


National Aeronautics and Space Administration 
Electronics Research Center 

Cambridge, Massachusetts, November 1968 
125-19-01-11-25 
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APPENDIX A 


PROGRAM LISTING 


1 1/13768 Bo /'BCfTTST PTO^'NOV- 0 OO 0 Ot 

♦JOB RSHIRLEY 

■ ' * 0 AT E-Tt NOVE MBER . 19-6- 8 

♦TITLE PROGRAM TO MEASURE THE DESCRIBING FUNCTION USING THE FFT 

♦A S S~IGN “ T = MT~2"A v2 =M T '1 8"»" 3=MT 1 A- 

♦ASSIGN 5 = CR 1 A * 6=LP 1 A 

♦•A-SS+G+T-^=C-P±A 

♦FORTRAN S * GO 

■e 

c 

C “TIA-fN "PROGRAM 


C 

C 

c 

•e 

c 

c 

c 

c 

c 

tr 



SENSE SWITCH FORMAT 

DO NOT ' SET "S ENSE S WTTC W ES DU R ING C OM P I LA t T O N - 

SET=SW I TCH IN. LIGHT ON 

IrOT-SET-SKfl-TCH- OOTi-Ll-GHT ' OF/ 


SWITCH NOV SET' NO/SET 

*************■*********#»*****♦*****#*♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦*♦♦*♦♦♦♦*♦♦♦♦*♦*♦*♦♦ 

1— * DO-N O r ~R € -A tr - &AT A— C A R DS ♦ "READ DAT A ' CA R DS-ANp— CALC a — TNRUT - 

2 * DO NOT CALC. INPUT * CALC. INPUT 

3-s. DO NOT PROCESS" DAT A ♦-PROCES-SHJATA 

4 ♦ DO NOT USE UNIT 3 TO SAVE DATA * USE UNIT 3 TO SAVE DATA 

5 * DO A DAT A— DUMP" ♦- L ND"DATA“DUMP“- " 

6 * DO PUNCH OUT DATA * DO NOT PUNCH OUT DATA 


- ***««■»**+«■* — HUMAN * *** ** *** * *** * SYSTE-M 

INPUT ERROR ♦ ♦ OUTPUT ♦ ♦ OUTPUT 

( PUT ) f ** * TEi » HUMAN * 10 N C ) » ■ VE H ICLE » 1 T W O T 

«•******■**# *■*■***■»***»* ************ *»************ 

*»* --■» OPERATOR — * DYNAMICS * ~ — -* 

* - * * * * * 

- * ******-******-** **-*#*■****■**■■**■* * 

* # 


C * *■ 

c ************************** ******* *********************** 

C*- — - 

c 

C" FORMAT ■ 'STATEMENTS 

100 FORMAT ( 1H 1 * 3X . 1HK . 6X . 8HFREO. ( < ) . 7X . 5HAE ( X ) . 10X .5HBE t K ) . 1 OX » 

rsHwrnn . iox.shbi m .i o x.sha 2 no .i o x .she tii o »a/. < is.teis.a) > 

101 FORMAT (///.4X.1HX.6X.8HFREQ. (X.) , 7X , 8HAR EMI ( X ) . 7X » 8HBREM1 (X) ,7X. 

1 8HAREM-2"(-K'7 v7X-,-8H&RE+42T-K-Tf A/-.-M-8-.-SE-15v4-n 

102 FORMAT ( 1E15. 8 ) 

10S-FORMAT-(-/r-/v//-.-4'X-,-THXr?<rX-,-8HFRED.-1-X-)-.-TX-,-6-HAR-tf-X:)-.7X-rTHPHAS-'t-Xr7-rfrX-? 

16HAR2 (X) .9X.7HPHA2 (X) .//.(IS.5E15.4)) 

104 FOR M AT ( 1 H 1.4X .1 H X.10X.4 H E (X ) » 9X ♦ 6HON C ( 10 > 9 X rO I I T W O ( X ) » / / t ( i 5 > 3 Cl D » 
24) ) 

•]:05-E©Rt4AT--t-H-Hr»4-X-»-lHXv9X-r6-HRU-T-(-X-Tr^X-.AHeH-X7-ri-irX-r6'HeN€-I-X-)-r^X"r6HTWO-(-k:-) 

1.//. ( I5.4E15.4) ) 

‘ 106 " FORMAT — ITH/r AXi~tHK‘."9X-r6'HPt>T-TX"T. Y-/"T"(""iE>TrlE _ l'3-* 4 Yj r 

107 FORMAT ( / / / . 1 OX , 8HPUTSQ = .1E15.4./.10X. 8HERRSO = .1E15.4.10X, 

114 H ER R 3 Q 2 PUTSQ - . 1 E 10 . 4 .A t IOX > 8 H ONCSO — .lC ^ 0 . 4 r 10 X ♦ 14 I I O N ESQ/ P UT 

2SQ = . 1E15.4./ .10X.8HTWOSQ = . 1 E 1 5 .4 , 1 OX . 14HT WOSQ/PUTSQ = .1E15.4) 
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n T 1 3 76 8 (TO/ 80 1.1 SI P ACE WO. 000002 - 

108 FORMATI///.4X.lHX,6X.8HFREO.<X>VfX»7HR£Ml<K)*iX.7HREM2<IC>*//» 

n r5T?FT5T¥n 

109 FORMAT (10X.3U0) 

T1TT F UHH A T~~r7~/Y il »XV5KJATirT V A X * 5 H J 8 I K > »4X «9BJC 1 X ) » / / . ( 2X i 1 9 ,4X » r3 v 4X T r3) 

1) 

- 111 FORMA T -tim-«rf 4X > lHK> 3 XVAH »tK r»//»TI9tE l» iAn 

112 FORMAT ( ///»4X» 23HREMNANT POWER AT ONE - .1E15.4./.4X, 

1 L iHKEWNAN I P UWEK "A I IWI3"~ »~71TI7'.'4 f 

113 FORMAT (///.4X,1HX.6X,8HFRE0.<K),7X.8HREMA1(K)»7X»8HREMA2(K) »//» 

1Tr5 - lfET5 -- 4yr — 

114 FORMAT <///.10X,$THE FOLLOWING REMNANT VALUES ARE CORRECTED FOR TH 

lE LOOP^GArNSV/Tf TV EVr7 , Kr - "ACT(JAL"i“PHI-WFA5r.~TTME^“ 1 + YPYC 'SQUARE 

2DS ) 

lib FORMA T ffeblO.4) 

116 FORMAT ( / »4X »8HCONS2 « .1E15.4) 

riT'FORMAT" 'r//V4Xrt"7HXT5TJ' = " VIFT5-.“4TllJX“.‘7HXTSUi . 1ET5V41 

C 

C-T51WE7*5T0W~5TXTEWENT3 

COMMON PUT<1080) ,E< 1080 > .ONE< 1080) .TWO (1080) »XONX. LOP. XANX. 

IN , H ;W.TH 1 NKTFUT5Q , hRR ' STT. ' ONESTl . 1 WOS ' U 

COMMON JA ( 29 ) . JB ( 29 ) . JC ( 29 ) ,WS < 4300 ) .XL < 29 > »KLE(29) 

tTIMEWSt 0 W'rC3TJlVAE:i X 5 T iHET I 5-r.^lTl‘5T-.'BTT T 5 T iX2 TI b T -.-B2 T l-5-TT 

1AREM1 ( 15) ,BREM1( 15) ,AREM2( 15) .BREM2I15) 

- - 0 - j MENSI on ~AR lT 1 5' ) .AR2TI 5 ) . PHA 1 <1 5 ) , PHA2U 5 >'7RERr( 1 5 ) » REM 2TI 5 ) 

DIMENSION AM (15) 

DIMENSION REmA1(15) ,REmA2<15) ' 

C 

c 

c THE PROGRAM CAN BE RECALLED TO THIS POINT BY IFINITIA AT ANY TIME 

C B y'hTt T ING f NT ER rU PT 33~A NO* TYPING" A~TArRIaGEWETUKN 

CALL fFlNITIA 

z 

C CONNECT THE CONSOLE 

S EOM'03li2'd* 

C 

C PUT THE AN A LOG C OMPU TE R I NTTF~YC MOOT - "' ' ' - - - 

CALL IC 


C SENSE SWITCH4 DETERMINES IF REWIND 3 OR NOT 

I F"< SENSE"SW" I TCH4f~ 98. 99 *' ~ -- - - 

99 REWIND 3 

98 CONTINUE " ' 

C 

C SENSE SWITCH1 DETERMINES IF FREQUENCY AND REGISTER VALUE CARDS ARE 

C _R_EAD _ _ __ 

IF f SENSE SWITCH!) 23.24 


C SENSE SWITCH2 DETERMINES IF INPUT CALCULATED 

23 IF f SENSE SWITCH2) 17.26 

C 

C READ IN 29 FREDuENCIES IN ORDER, HIGH FREQUENCIES FIRST. 

r TN''RADTANS7srCOND',''STARTmG"liTTH‘A‘‘RtHNANrTTEEaUENCY _ 1RPJD _ 'ALTERNA- 

C TING THEREAFTER WITH THE INPUT FREQUENCIES. THE FREQUENCIES MUST 

Z ATU'B'E''rNTFGER"'W0i:TTPrET"UF“2TPrr7R0NniiifF‘.‘'AND I HE""NUMBER 'OF 

C SECONDS PER CYCLE MUST EQUAL 41 (DELTA) » WHERE I IS AN INTEGER AND 

Z DELTA IS Th£ - T IME BETWEEN DATA P0INT5. 

24 READ! 5 .102 ) < W < X > »X-1 .29 ) 
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11/13/68 80/8U LIST PAGE NO. 

WR I TH ( 6.111) (K.W(K) ,K»1,29) 

C 

C READ IN THE RtGISTER VALUES FOK T Ht AUUiTIVt PART OF THE FFT 

READ (5«109) ( JA(K) ,JB(K) ,UC(K > .K-1,29) 

WRITE! 6. 110) < JAU) ,JU(K) ,JC(K) ,K»1,29> 


INPUT (CALCULATED EVERY DELT) 

CONS1 SCALES THE INPUT 
26 C0NS1=9.5 

REWIND INPUT TAPE PRIOR TO STORING INPUT 
REWIND 1 

ONW SETS THE STARTING TIME FOR THE INPUT. AND IS NEGATIVE 
SO THAT THE HUMAN OPERATOR IS IN A STEADV-STATt TRACKING 
CONDITION WHEN THE ONSET OF DATA-TAKING OCCURS 
ONW=-1080. 


DELT IS THE TIME INCREMENT BETwEEN INPUT VALUES. AND MUST EQUAL 
DELTA, THE TIME INCREMENT BETWEEN DATA POINTS. 


DELT= • 02 

THE AM ( K ) SCALE THE INPUT SINUSOIDS 

AM ( 1 ) = .2 

AM ( 2 ) =-• 2 

AM ( 3 ) = • 2 

AMI 4) =-.2 

AM ( 5 ) = « 2 

AM ( 6 ) =-l , 

AM ( 7 ) = 1 . 

AM(b)=-l. 

AM ( 9 ) = 1 . 

AM<10)=-1. 

AM ( 1 1 ) = 1 . 

AM( 12) =-l. 

AM ( 1 3 ) = 1 • 

AMI 14) =-l. 

AM ( 1 5 ) = 1 • 


C 

C 

C 


the input is calculated in blocks of 540 values and stored on 
magnetic tape (unit i>. 

DU 19 J= 1 » 30 
UO 1 K = i *540 
PUT(K) =0. 

ONW=ONW+l • 


T =ONW*UE LT 
DO 20 L= 1 . 14 

PUT ( K ) =PUT ( K ) +AM ( L ) *S IN ( T*W ( 2*1 ) ) 

20 CONTINUE 

PUT ( K ) =CQNS1*PUT ( K 1 
1 CONTINUE 
M= 1 
N = 540 

CALL BUFFERUUTI 1 ,1,PUT(M) ,2*(N~M+1) , I STATUS) 
CALL GOTO! ISTATUS) 

19 CONTINUE 


17 CONTINUE 

THE INTERRUPT I S. CONNECTED » BUT NOT ENABLED 


OOOOU3 
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11/13/68 80/80 LIST PAGE NO. 000004 

CONNECT (40.INTR) 

EOM 020020 
POT =00700000 


INITIALIZE FOR THt RUN 

ZERO THE BUFFER AREAS 
DO 92 1.1080 

PUT(K)=0. 

E ( K ) = 0 - 
ONE ( K ) =0. 

T WO ( K ) =0 • 

92 CONTINUE 

ZcRO T He DATA TAPE 
REWIND 2 
DO 91 K = 1 » 90 
N = 640 
M= 1 

CALL BUFFEROUT (2,l,t<M),2*< N-M+l ). I STATUS ) 

CALL GOTO! ISTATUS) 

91 CONTINUE 

THE MAGNETIC TAPE UNITS ARE INITIALIZED, UNIT 1 FOR THE INPUT, 
AND UNIT 2 FOR THE DATA. 

REWIND 1 
REWIND 2 

ZtRO THt RtGlSTERS WhtRt THt FOURIER COEFF I t I tNTS ARE TO 
BE CALCULATto. 

DO 31 K = 1 , 1 5 
At ( K ) = 0. 

BE ( K ) = 0. 

A1 ( K ) =0. 
b 1 ( K ) = 0-. 

AZ ( K ) = 0. 

□ 2 ( K.) =0. 

ARtMl ( K.) =0. 
dRtMl ( K) =0. 

ARtM2( K ) =0. 

BREM2 ( K ) =0. 

31 CONTINUt 

INITIALLIZt FOR Ti-it ADDITIVt PART OF THt F F T 
DO 10 NOW= 1,29 
K.L 1 NOW ) =-l 
10 CONTINUt 

DO 14 NOW= 1,4300 
WS ( NOW ) = 0 . 

14 CONTINUE 

INITIALLIZt THE CuuNTERS FOR THt RUN 

K.R I N K. IS A COUNTER TO Dt TERM 1 Nt THE LOCATION FROM WHICH THE 
NtXT INPUT VALUE SHOULD bt TAKEN FROM 
KR I NK= 0 

KONK COUNTS THE INTERRUPTS, DETERMINES WHEN Tnt ONSET OF DATA- 
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C TAKING SHOULD OCCUR* ANU WHEN OATA-TAKING IS COMPLETED. 

KONK = 0 

C KANK IS THE HALF REGISTER COUNTER* 1 TO 5A0 

KANK=0 

C MO DETERMINES WHICH HALF 00 THe INPUT ooFFlR IS BeING USED 

M0 = 0 

C 8 IS A COUNT eR UN THE DATA USED DURING THE DA T A-PROCES S I NO 

b = -l . 

C LOP IS THE Flag SET BY INTR TO END DATA-TAKING 

LOP = 0 

C PUTSQ AND ERRSO ARE THE INTEGRAL SQUARE INPUT AND ERROR 

put sq= o. 
errsq=o. 

ONESQ=0. 

TWOSQ = C. 

xisu=o . 

X2SU=3. 

INITIALLIZe INPUT bUFFtiR FOR TnE RUN* I.c., FILL uOTn HALVES 
WITH INPUT VALUES 
M= 1 
N = b A 0 

CALL BuFFeRINI 1* 1 , PUT ( M ) ,2*<N-H+1 ) , ISTaTuS) 

CALL GOTO! ISTAT'JS) 

M = 5 A 1 
N = 1 0 80 

CALL BuFFeRINI 1*1*PUT(M)*2*( N-i'i + i ! > I ST ATOo ) 

CALL GOTO ( ISTATUS) 

WAIT TO START RUN ON SIGNAL FROM The OPERATOR 
PAUSE 


511 CONTINUE 
SKS G 3 0000 
oRU 5 1 3 S 
GO TO 511 
513 CONTINUE 

PUT TnE ANALOG COMPUTER INTO COI-’PuTe mode 
CALL COMPUTE 

enable the interrupt 
EOM 031032 


11 CONTINUE 

C CHECK TO SEE IF IT IS THE END OF DATA-TAKING 

IF (LOP.EO.l) GO TO 2 

C IT IS NOT THE END OF DATA TAKING. WAIT FOR INTERRUPT 

GO TO 11 
C 
C 

C IT IS TnE lND OF DATA TAKINO* uu ON 

2 CONTINUE 

C TURN OFF THE INTERRUPT 

S EOM 031033 

C PUT THE ANALOG COMPUTER INTO TH t HOLD MODE 


OCGOGd 
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CALL HOLD 

C TAKE THE INTEGRAL SQUARE MEASURES 

CALL ADL ( A,PUTSU»ERRSU,0NESQ,TW0SQ,X1SU,X2SQ> 

C PUT THt ANALOG COMPUTtR INTO IL MODt 

CALL 1C 


SENSE SWITCH3 DETERMINES WHETHER TO PROCESS THE DATA AND TYPE 
THE RESULTS. OR WHETHER TO RE- I N I T I ALL I ZE FOR THE NEXT RUN 
IF (SENSE SWITCH3) 17,22 

' DATA PROCESSING IS DESIRED, GO ON 
22 CONTINUE 

THE FOLLOWING SAVES THE LAST 5*.0 DATA POINTS 
N= 1080 
M = 5 A 1 

/ CALL BUFFERQUT ( 2 ,1 »E(M) ,2*<N-M+1 ) , JSTA I US ) 

CALL GOTO ( ISTATUS) 

CALL BUFFEROUT(2,l,ONE(M) ,2* (N-M+l) .ISTATUS) 

CALL GOTO( ISTATUS) 

CALL BUFFEROUT ( 2 ,1 ,TWO(M) ,2* (N-M+l ), ISTATUS) 

CALL GOTO ( ISTATUS) 


THE FOLLOWING PERMITS A TOTAL OR PARTIAL TAPE DUMP 
IF (SENSE SWITCH5) 89,90 
89 REWIND 1 
REWIND 2 
N = 540 
M = 1 

CALL BUFFER I N(.l ,1 ,PUT|M) ,2* (N-M + l ) , ISTATUS) 

CALL GOTO ( ISTATUS) 

N= 1080 
M = 5 A 1 

CALL BUFFERIN ( 1,1 .PUT (M) ,2* l N-M + l ) .ISTATUS) 

CALL GOTO ( I STATUS ) 

WR I T E ( 6 . 106 ) (J,PUT (J) ,J=1,1080) 

C GO PAST THE INITIAL SPURIOUS DATA POINTS (CAUSED bY INTERRUPT 

C ROUTINt). 

M = 1 
N = SA0 

DO 16 J=l,3 

CALL BUFFERIN ( 2, 1 ,E (M) ,2# (N-M + l I, I STATUS I 
CALL GOTO ( ISTATUS) 

16 CONTINUE 

C NO DETERMINES EXTENT OF THE DUMP 

NO =3 

DC 88 K = 1 , NO 

CALL BUFFERIN ( 1,1 »PUT(M> ,2* (N-M+l ) .ISTATUS) 

CALL GOTO ( ISTATUS) 

CALL BUFFER IN ( 2, 1 »EM) ,2* (N-M+l ), ISTATUS) 

CALL GOTO ( ISTATUS) 

CALL BUFFERIN! 2, 1 ,ONt(M) ,2*(N-.i+l ) , ISTATUS) 

CALL GOTO ( ISTATUS) 

CALL BUFFERIN ( 2, l.TWO(M) , 2* ( N-M + l ) , 1ST ATU8 ) 

CALL GOTO ( ISTATUS) 

L= ( K-l )*5A0 
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WRITEI6.105) ( IL+Jl »PUT( J).E( J) »OHE( J1 »TWO( J) »J«1.540> 

88 CONTINUE 

90 CONTINUE 

C 

C 

C~~ " PROCESS' THE DATA FOR 10.800 POINTS 
C 

c 

C MULTIPLICATIVE PART OF THE FFT USED TO CALCULATE THE FOURIER 

C COEFF ICl ENTS 

C __ _ _ _ 

C CALCULATE 'tHE' FOUR IER COEFFICIENTS F0R~THE~1 ?EWIAnT 
DEZ*3. 1415927 

DO 21 K-'1*29.Z” ’ 

NEMK.+ 1) /2 

U£BT ij B T RT +J A T RT-T 

LEE1»JB(io+2*JA<IO-1 

rEB2 VJB ( kT+2*'JATRT-T "" “ " 

LEE2* JB ( K)+3*JA( K) — 1 

lEM1=lEE1 _j4 (k j /2 

LEM2=LEE2-JA(K>/2 

a^remITneiVoV 

AREM2 ( NE )*0. 

B'REmITNET“6V 

BREM2 ( NE ) =0. 

jst _ ( -j A<k , / 2 )-\ 

00 27 JI =1 * JST 

*UOU»FLOAtTjI ) 

VOU=FLOAT ( JA ( K . ) ) 

T OU* UO 07 VOU' "• 

SINUS«S^N(DE2*T0U) ' 

C0SUS=S(jRm.-SINUS**2. J = 

AREM1 < NE J *ARcMytAt£ Jj+< WSU.EB1+ JI >>WS( LEEl-JI ) »*SI MS 
_____ ( __ _- REf(|2 ( P4£ j ( LEB2+ JJ t^-WSl LEE'2- j I ) ) *S I'NUS' 

BREM1(NE)=BREM1(NE) LEE 1-JI ) )*COSUS 

0-R- E ^2TNEr=BREM2'(N'E" ) + (WS(LEB2+jf)-WS(LEE2-Jl D*C<5SUS'' 

27 CONTINUE 

AREM1(NE)=AREM1(NE H*IS<LEM1 ) 

AREM2 ( NE ) = AREM2 < NE )+*S( LEM2 > 

S R EM IT NE Y=~ B REM ITNE ) CE ET 1 

SREM2 ( NE ) =BREM2 ( NE ) -WS< LEE2 ) 

2TT0NTTNUE: 

C 

Z Calculate The FOURIER COEFFICIENTS FOR the DESCRIBING FUNCTION 

00 28 K=2 » 28.2 

NE=K./~ 2 “ * ' ‘ 

LEBE»JB()0-1 

LEEE* JB fk) +JA t K.T-1' ' ' ' - — - ' 

LEBl*JBm+JA<K)-l 

LEE1*JB<K)+2*jA(K>-1 

LEB2sJB( (0+2*JA(<)-l 

L EE2 = JF( K > +' 3*~JA ( kT-I ' 

L£ME=*LEEE-JA(K.) /2 

L E M 1 = LEEl'- J A ( kT7 2 " " " ' 

LEM2=LEE2-JA(K)/2 
AE < NE ) *0. 

BE(NE) -0. 
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DO”32. «C»l,'l4 

DENOM“AE I K )^AE ( K ) +BE ( K ) *BE ( K ) 


R1 OUTSORT! (A1(IC)*A1(IU+&1(IC)*B1(IC) >/l 
AR2 ( K ) “SORT ( IA2 ( K) *A2 ( ICl +B2 ( IC)*B2 ( (C ) ) /DENOM1 


PHA1 ( tC ) »57« 3* ( ATAN2(B1(K) » A1 ( K. ) )- ATAN2 I BE ( 1C) » AE ( K) ) ) 

PHA2(K) = 57.3*< ATAN2(B2<K).A2(K.U-ATAN2(BE(IC> ,AE(IO ) ) 

USE THE ASSUMPTION THAT THE HUMANS PHASE LEAD IS "LESS THAN" 
180 DEGRESS TO CORRECT FOR THE LOSS OF PHASE INFORMATION 


IN THE ARC-TANGENT ROUTINES 
IF (PHAl(K).LT.ieO.i GO TO 93 
"" P'HAl'l K") = PHA1 <K) -3F6V ' “ 

93 IF < PH A2 (IC).LT.180.) GO TO 94 

P'H A 21 X 1 "* PHAS ( <T - 36 0 • " 

94 CONTINUE 


32 CONTINU 


CONS2 SCALES THE REMNANT 
CONS2* ( DELT) **2/4( P I ) T 
DELT = .02»"~T*216 
CONS2*1.47E-7 


REM1 ( K ) *CONS2* ( AREM1 ( K ) *AREM1 C K.I+BREM1 1 K ) *BREM1 C 1C) ) 

ITEW2TKT*C‘0NS2*T7fREFr7TTCTTFATTEH2TXT+BREH2TiCT*FREM2TKTy 

33 CONTINUE 

TflTER'F0CATF"F0R""TFF"RFfWA'N7"~AT'THFTNPUr“nTEUUENC I ES 

DO 230 IC«1.14 


( K I =REM1 ( 1C) + ( REM1 ( K+l)— REM1 1 1C) )*()»( 2*IC ) — W l 2*IC— 1 ) )/( 
1W( 2*IC-1) ) 
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REMA 2 ! K) =RSM 2 ( K)-M REM 2 ( K+l )-R£M 2 1 IQ ) » ( W ( 2 «K ) -W ! 2 »K -1 ) ) / 1 W ( 2 »K +1 ) - 

' 1W<2#K-1 )') 

230 CONTINUE 

C CALCULATE THE REMNANT "POWER ' ' 

REMPW 1 = 0 • 

REMPW 2 = 0 . 

DO 231 K=l ,14 

R'EMPW 1 =REMPW 1 + REMA 1 <K)*(W< 2 *K- 1 )-W< 2 *K + 1 > ) 

REMPW2 = REMPW2+REMA2 (X )* ( W ( 2*K- 1 ) -W ( 2*K+1 ) ) 

231 CONTINUE 


WRITE OUT THE FOURIER COEFFICIENTS OF THE SYSTEM ERROR (E). THE 
HUMANS OUTPUT (ONE), AND OF THE SYSTEM OUTPUT (TWO). 

WRITE! 6.100) (K,Wl 2*K.) . AE ( K ) » BE ( K ) , A1 ( 1C) >B1 ( K. ) . A2 ( < ) . B2 l K ) »K.= 1 . 1 A ) 

WRITE OUT THE FOURIER COEFFICIENTS OF THE REMNANT AT THE HUMANS 
OUTPUT (ONE). AND AT THE SYSTEM OUTPUT (TWO). 

WRITE l 6. 101) (K.W(2*K-1> » AREM1 ! K ) » BREM1 ( K ) » AREM2 l K. ) ,BREM2(K> »H=1.1 

15) 

WRITE OUT THE HUMAN OPERATORS DESCRIBING FUNCTION AND REMNANT AS" - ' 
WELL AS THE SYSTEM OPEN LOOP DESCRIBING FUNCTION AND REMNANT 
WRITE(6»108) ( K»W( 2*K-1 ) . REM1 ( K ) .REM 2 ( K j »K=1 , 15 ) 

WRITE(6.116) CONS2 
‘ C WRITE OUT THE' HUMANS REMNANT' AT THE' INPUT FREO'DE'N'CITS 

WRITE (6,113) (<»W(2*K.) .REMAl l<) »REMA2t)0 »K = 1 » 14 ) 

C 

WR I TE ( 6,103) <K,W< 2*K) ,AR1(K) .PHAl(K) ,AR2(K.) ,PHA2U) »<=1 .14) 

DO 239 K=1.14 

REMAl ( O =REMA1 ( K.)* t (l.+AR2(K) )**2.) 

REMA21 <r=REMA? ( k) »('( 1.+AR2 ( K ) ) **2 • ) * ( AR1 ( K. 1**2. > /TAR2 t K‘) #*2 . ) ~~ 

239 CONTINUE 

WR ITE( 6 , 1 1 A ) 

WRITE! 6,113) (K,W(2*K) .REMAl (K) »REMA2(K)»K=1»14) 

WR ITE( 6,116) CONS2 
C WRITE OUT THE ERROR SCORES 

errsy=er r s q / p u f s Q “ ' " 

ONESY=ONESQ/PUTSO 
T WOSY = TWQSQ/PUTSQ 

WRITE (6, 107) (PUTSQ.ERRSQ.ERRSY ,ONESQ .ONE5Y » TWOSQ , TWOSY ) 

WR I TE ( 6 » 1 17 ) X1SQ.X2SQ 
PUTSO=PUTSQ/216. 

ERRSQ=ERRSQ/216, 

ONESQ=ONESQ/2 1 6. 

TWOSQ= TWOSQ/ 216. 

WRITE! 6,107) ( PUTSQ ,ERRSQ .ERR SY .ONE SO .ONESY , TWOSQ » TWOSY ) 

C WRITE OUT THE REMNANT POWER 

WR I TE ( 6, 1 12 ) ( RE MPW1 .REMPW2 ) 

REMPW1 =0 • ~ 

REMPW2 =0. 

DO 240 K.= l ,14 

REMPW1=REMPW1+REMA1 (K)#(W( 2*K-1 )-W( 2#K + 1 ) ) 

REMPW2»REMPW2+REMA2(K)*(W('2*K-l)-W(2#k+i) ) " 

240 CONTINUE 

WRITE! 6. 112) ( REMPW1 .REMPW2 ) ' 

C SENSE SWITCH 6 DETERMINES WHETHER TO PUNCH OUT THE ANSWERS ON CARDS 


23 


nonnrinn r\ n n n n n n n n 


PAGE NO. 000010 


11/13/68 80/80 LIST 

IF (SENSE SWITCH6) 35,36 
PUNCH OUT THE DATA ON CARDS* 14 CARDS WITH AR 1 ,PHA1 , AR2 , PHA2 . 
REMA1, AND REMA2 * PLUS 1 CARD WITH PUTSU* ERRSO.ONESU * TwOSCl » 
REMPW1, AND REMPW2. 

35 WRITEf 7,115) ( AR 1 ( K) , PHA1 ( K ) ♦ AR2 < K.) , PHA2 ( X) *REMA1 ( K ) ,REMA2( K.) * X= 1 , 
114) 

WRITE! 7,115) ( PUTSQ , ERRSG »ON£SQ » TWOSQ *REMPW 1 ,REMPW2 ) 

36 CONTINUE 


THE FOLLOWING STORES EITHER ONE (M) OR TWO(M) ON UNIT 3 FOR 
LATER PROCESSING FOR THE REMNANT 
SENSE SWITCH4 DETERMINES WHETHER TO SAVE ONE ( M ) OR TWO(M) 
IF (SENSE SWITCH4) 95,97 
97 REWIND 2 
REWIND 1 


M= 1 
N = 540 

THE FOLLOWING AVOIDS THE SPURIOUS DATA POINTS 
DO 521 K=1 ,3 

CALL BUFFER IN ( 2,1, TWO(M) ,2*(N-M+1 ) , I STATUS) 
CALL GOTO ( 1 STATUS) 

521 CONTINUE 

CALL BUFFERIN(1»1,PUT(M) ,2*<N-M+1 ) , I STATUS) 
CALL GOTO ( ISTATUS) 

CALL BUFFERIN( 1,1,PUT(M) ,2*(N-M+1 ) .ISTATUS) 
CALL GOTO ( ISTATUS) 

DO 96 K= 1 , 20 

CALL SUFFER IN ( 1,1,PUT(M) ,2* (N-i'i+1 ) .ISTATUS) 
CALL GOTO! ISTATUS) 

CALL BUFFER IN ( 2,1,E(M) ,2*(N-M+1 ) .ISTATUS) 

CALL GOTO ( ISTATUS) 

CALL BUFFER I N ( 2,1, ONE ( M ) ,2*(N~M+1 ) , ISTATUS) 
CALL GOTO ( ISTATUS) 

CALL BUFFZRINC 2 , 1 , T WO ( M ) ,2#(N-M+1 ) , ISTATUS) 
CALL GOTO! ISTATUS) 

CALL BUFFEROUT (3,1, PUT ( M ) ,2*(N-M+1) .ISTATUS) 
CALL GOTO ( ISTATUS) 

CALL BUFFEROUT ( 3,1,E(M) ,2*(N-M+1) , ISTAlUS) 
CALL GOTO ( ISTATUS) 

CALL BUFFEROUT ( 3,1 .ONE (M) ,2*<N-M+1) .ISTATUS) 
CALL GO TO ( ISTATUS) 

CALL BUFFEROUT ( 3,1 , TWO! M) ,2*(N-M+1> .ISTATUS) 
CALL GOTOt ISTATUS) 

96 CONTINUE 
95 CONTINUE 


RETURN TO INITIALIZE FOR THE NEXT RUN 
GO TO 17 
200 STOP 


INTERRUPT SUBROUTINE (INTERNAL) 
I NTR SERVICES THE INTERRUPT 
SUBROUTINE INTR 
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KRINK IS THE REGISTER COUNTER. 1 TO 1080 
KR I NK=<R I NK+ 1 

KONK. IS THE TOTAL COUNTER . 1 ON UP 
KONK=<ONK+l 

RANK. IS THE HALF REGISTER COUNTER. 1 TO 500 
KANK = K*N.<+1 . 

IF (KONK. GT. 1060) GO TO 1 
'l IF (.<ANK.GE.541 ) GO TO 3 
4 CALL DAL(O.PUT(<RINK) ) 

12 RETURN 

1 } F (KONK.GT. 11880) GO T'o 5 
IF ( <A^K.pE.541 ) GO TO 6 
8 CALL DALIO.PUT (KRINK) ) 

CALL ACL (O.E( KRINK) ,0n£( KRINK) . T w.O ( KR I NK > ) 

ADDITIVE PART UF THE FAST FUUKItR TRANSFORM 

KL = A COUNTER. J TO ( 2* J A ( <) * JC ( < ) - 1 > 

JA=HALF THE NUMBER OF CYCLES PER SECOND/i) i V I SOR 
Jd=REG I STER START FOR EACH FREQUENCY 

JC = D I V I SOR USED ON HALF THE NUMBER OF SAMPLES PER CYCLE TO GET JA 

ADoITIVE PART OF THE FFT 

ADD Trie DATA INTO ThE WS(.<) 

00 27 <=1,29 

23 KL ( K ) =KL ( < ) + 1 
JNK=JA (<)*JC( K> 

IF (KL I.O.GE.JNK) GO TO 32 
<LE ( .<) = KL ( <) / JC ( < ) + Ja(K) 

Go TO 24 

32 Jm,< = 2#JA(,<)*JC(<) 

IF(<!_<<) .6E. JAK) GO TO 33 

<LE (<> = ( <L (.<)/ JC (< ) ) - JA ( K ) + Jb ( K ) 

GO TO 25 

33 K. L (.<)=- 1 
GC TO 23 

24 Ju=KLE ( < ) 

Jc=<LE ( < )+JA( < ) 

JF = KLc. ( < ) + 2* JA ( .< ) 

wS ( JC ) =WS ( JC ) + E ( <R 1 1\< ) 
wb ( JE ) = «S ( JE 1 +ONE ( <R I NK ) 

*S< JF ) =WS ( JF ) + T WO (KRINK ) 

GO TO 26 

25 JD=KLE(K) 

Jt = <LE <K)+JA( <) 

J.-=KLE ( K ) + 2 * J A ( .<) 

WS ( JO ) =WS ( JD ) -E ( <R I NK ) 
wS < JE ) =WS ( JE ) -ONE ( KR I NK i 
as ( JF ) = WS ( JF )~TWO( KRINK) 

2o continue 

27 CONTINUE 
RETURN 
3 KAN<=: 

IF (MO.GE.l) GO To 7 
MO= i 
N = 540 
M=1 


00O011 
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9 CALL BUFFERIN( 1,1. PUT(M) ,2»(N-M+1 ), ISTATUS) 

GO TO A 

7 M0=0 _ _ _ _ - 

N=1080 
M = 541 

K.R INK= 1 ' 

__ GO TO 9 

6 KANK=1 

IF (MO.GE.l) GO TO 10 
M0= 1 
N = 540 
M = 1 

11 CALL BUFFER INt l.l.P UT(M) . 2* ( N-M +l ) , 1ST ATUS ) . 

CALL BUFFEROUT ( 2 ,1 ,E (M) *2* ( N— M+l ) * 1ST A t US j 
CALL BUFFEROUT (2.1. ONE (M) »2*(N-M+1> .ISTATUS) 

CALL BUFFEROUT (2.1.fW0(M) . 2* ( N-M+ 1) .ISTATUS ) 

GO TO 8 
10 M0=0 

N = 1080 _ 

M = 641 

KR I NK= 1 
GO TO 11 
5 LOP= 1 
GO TO 12 


SUBROUTINE GOTO (ISTATUS) 

c Subroutine (internal) to handle tape read and writes 

7 GO TO (6,4,5,5.51 ISTATUS 
6 GO TO 7 

5 WRITE (102,200) ISTATUS 

4 RETURN : 

200 FORMAT ( SBUFFER I N STATUS WORD =S , 1 2 ) 

END 

*LOAD X 

*DATA 

39.269908 

26.179938 

19.634954 

15.707963 

13.089969 

8.7266460 

7.8539816 

6.5449847 

5.2359877 

4.3633230 

3.1415927 

2.6179939 

2.1816615 

1 . 7453292 

1 . 5707963 

1.3089970 

1.0471975 

. 8 7266462 “ 

. 7272 2052 

. 58177642 ' ^ 

.52359877 
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.43633231 




.34906585 

.29088821 

.26179938 

.17453292 

. 14544410 
.11635528 

.08726646 

4 

1 

1 


6 

13 

1 


8 

31 

1 


10 

55 

1 


12 

85 

i 


16 

121 

i 


20 

175 

l 


24 

235 

l 


30 

307 

l 


36 

397 

i 


50 

505 

i 


60 

655 

i 


72 

835 

l 


90 

1051 

l 


50 

1321 

2 


50 

1651 

3 


60 

1471 

2 


60 

rsoi 

3 


54 

1981 

4 


54 

2143 

5 


50 

2305 

6 


60 

50 

2455T 

2635 

9 


60 

TTS5 

9 


60 

2965 

10 


108 

3415 

10 


90 

3145 

M _ 



3739 

13 ~ 


90 

4009 

20 
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APPLICATION OF A MODIFIED FAST 
FOURIER TRANSFORM TO CALCULATE 
HUMAN OPERATOR DESCRIBING FUNCTION 

By Richard S. Shirley 
Electronics Research Center 


ABSTRACT 

A modified fast Fourier transform (FFT) is used in a hybrid 
computer program to permit processing of tracking data during a 
run to yield the human operator's describing function almost im- 
mediately after the data-taking period. The computer processing 
time is substantially reduced at no cost in accuracy. 
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